One-way quantum computing in optical lattices with many atom addressing 
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One of the fundamental conditions for one-way quantum computation (1WQC) is the ability 
to make sequential measurements on isolated qubits that comprise the highly entangled resource 
for 1WQC, the cluster state. This has been a significant impediment in the implementation of 
1WQC with ultracold atoms confined in optical lattices, because the width of the measuring lasers 
is generally much greater than the atomic (qubit) spacing. We demonstrate that deterministic 
1WQC is nevertheless possible, with a polynomial increase in the number of operations, as long as 
the center of the beams can be positioned with high accuracy. Extending the number of cluster 
atoms, the scheme is also able to compensate for accidental measurements of an arbitrary number 
of nearby qubits. 

PACS numbers: 03.67.Lx, 03.75.Ss, 05.30.Jp 
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I. INTRODUCTION 

In one- way quantum computing (1WQC), first pro- 
posed by Rausscndorf and Briegel [1| , the computation 
proceeds entirely by performing a sequence of single qubit 
measurements on a special many qubit entangled state, 
known as a cluster state. This approach is wholly equiv- 
alent to the quantum circuit model, but provides two 
distinct practical advantages. First, all of the required 
entanglement is generated at the outset in only two steps, 
involving all of the particles simultaneously. Second, no 
entangling gates are required between arbitrary (and po- 
tentially distant) qubits during the computation. The 
1WQC approach therefore lends itself to any physical ar- 
chitecture where one can induce genuine multipartite en- 
tanglement and perform singlc-qubit measurements. To 
date, 1WQC has been demonstrated in linear optical sys- 
tems by the implementation of Grover's search algorithm 
for up to four qubits 0]. 

Ultracold atoms confined in optical lattices arguably 
constitute the most promising test beds for quantum 
simulation and computation, and are natural candidates 
for the implementation of 1WQC. This optimism is due 
to the exceptional control of long-lived atomic internal 
states, the lattice parameters, and the interactions be- 
tween particles. Under the right conditions, exactly one 
bosonic atom will occupy each site of the optical lat- 
tice at low temperatures as the gas undergoes a quantum 
transition into a Mott insulating phase [3|. By adjusting 
the orientation and strength of the applied lasers, the 
Mott phase can be realized for atoms in effective one- 
dimensional (ID), 2D, or 3D lattices Because two 
internal states can be chosen as computational registers, 
each atom corresponds to a qubit. In the Mott phase, 
these (neutral) qubits are effectively stationary and non- 
interacting. The massively entangled cluster state that 
is the central resource for 1WQC can then be generated 
with relative ease @, H, 0| : starting with a single large 
quasi-2D array of atoms [1, Q, entanglement between 
neighboring atoms can be done in parallel with state- 



dependent collisions [Hj, or with tunable spin-spin 
interactions [H, HU. Allowing these operations to oc- 
cur for the right amount of time on qubits initialized 
in the zero eigenstate of the Pauli X operator, one can 
implement (up to local unitaries) the maximally entan- 
gling controllcd-phase (CZ) gate that generates the clus- 
ter state. 

While the relative ease with which one may gener- 
ate large cluster states containing thousands to millions 
of physical qubits appears to strongly favor ultracold 
atoms in optical lattices, these systems tend to suffer 
from one important shortcoming: the practical difficulty 
of measuring single qubits. Most optical lattice exper- 
iments to date utilize 87 Rb atoms (H, [ll|, with a reso- 
nant absorption at a wavelength of 780 nm. The wave- 
lengths of the lasers forming the optical lattice are gen- 
erally chosen to be sufficiently close to this resonance 
(usually A w 800 nm) in order to yield strong con- 
finement while minimizing spontaneous emission which 
causes heating and loss of atoms. Assuming perfectly 
counter-propagating lasers, the spacing between sites of 
the optical lattice is therefore a w 0.4 /mi. Yet the ad- 
ditional lasers that would be used to apply rotations or 
measurements on selected qubits can generally be focused 
to widths on the order of 3 — 4 /im, which is a few times 
the fundamental resolution limit. 

Much effort in recent years has been devoted to im- 
proving the addressability of single qubits in optical lat- 
tices. An experimental scheme has been realized wherein 
atoms occupying every site of an optical lattice are trans- 
ferred to every third site [Til ]; in principle this process 
can be repeated until the desired separation is reached. 
A CO2 laser generates an optical lattice with the long 
spacing a = 5.3 /im (l5j : while a Mott transition in this 
lattice is difficult to achieve, unit-filling in a finite re- 
gion may be reached by physically moving atoms through 
state-dependent laser manipulation fl6| |. This kind of 
state-dependent transport can also be used to move a 
small number of 'marker' atoms through the lattice, gen- 
erating entanglement with the cluster qubits (l7l. Ha, fl9| : 
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measurements of the (well-separated) marker atoms can 
then be unambiguously made. A recent theoretical pro- 
posal makes use of the energy shifts induced by additional 
lasers to ensure that only one of the atoms is resonant 
with the measuring laser (20| . Another uses interference 
patterns generated by multiple additional lasers oriented 
at different angles to localize atoms [2l[ ; unfortunately, it 
is generally not feasible experimentally to generate beams 
from very many directions or to have the attendant op- 
tical access. A very recent proposal is to pump atoms 
in the vicinity of the target into internal states that are 
decoupled from the subsequent manipulations [22| . 

In the present work we show that high-fidelity 1WQC 
does not in fact require single-qubit measurements, and 
in principle is experimentally feasible with lasers that si- 
multaneously impinge on a large number of atoms from 
the same direction. The spatial variation of the Gaus- 
sian measuring beam is the crucial ingredient making this 
possible. In general, projective measurements on multi- 
ple qubits generally yield mixed states because the lack 
of spatial resolution prevents knowing which qubits were 
projected into the computational basis states |0) and |1). 
In order to reduce the uncertainty, one could simply re- 
measure different subsets of the qubits. In this scenario 
one would expect to require TV — 1 additional measure- 
ments for each initial projection of TV qubits to reduce 
the uncertainty to zero and yield a pure state. An even 
simpler approach would be to arrange that the states of 
all qubits to be irradiated (except the one of interest) 
are in a simple fiducial state, such as one of the com- 
putational basis vectors, so that the ambiguity in the 
result of the multi-qubit measurement is always zero. In 
this second scenario, one would expect TV — 1 additional 
unitary operations in advance of each measurement. As 
discussed in detail below, it is always possible in principle 
to implement the second approach, ensuring high-fidelity 
1WQC even for large numbers of simultaneously mea- 
sured qubits. 

The manuscript is organized as follows. The basic 
ideas of 1WQC and the fundamentals of implementing 
rotations and measurements for atoms in optical lat- 
tices are reviewed in Section [TTJ The details of how 
to initialize atoms in the states required to form one- 
dimensional cluster states with many atom rotations are 
described in Section IIII1 and the protocol for perform- 
ing (non-universal) 1WQC with multiple measurements 
of atoms on a line is discussed in Section IIVI The full 
two-dimensional protocol is described in Section [V] The 
results are summarized in Section [VTl with comments on 
how the work could be extended. 



II. 1WQC AND SINGLE ATOM OPERATIONS 

A. One-way Quantum Computing 

Starting with a periodic two dimensional array of TV 
qubits, a cluster state can be prepared by first initializing 



all TV qubits to the +1 eigenstatc of the Pauli TV opera- 
tor |+) = (|0) + |l))/\/2, followed by the application of 
controlled-phase CZ = diag(l, 1, 1,-1) gates between all 
nearest neighbors. Computational measurements, made 
in the eigenbasis of the Pauli Z operator, are then per- 
formed to remove qubits from the initial cluster state 
and to enable the construction of a computation-specific 
cluster state (hereafter denoted a computational clus- 
ter state) with the desired connections. In the simplest 
(though not most computationally efficient) approach, 
the resulting graph state is a spatial representation of a 
quantum circuit: horizontal chains of entangled physical 
qubits encode a single computational qubit, and vertical 
links represent two-qubit entangling operations. 

Actual computation proceeds via quantum teleporta- 
tion by sequentially measuring the state of each qubit in 
a basis defined by the states |± € ) = (|0) ± e~^\l)) /V2. 
Alternatively, the state of each qubit is rotated by an 
angle £ around the Z axis, followed by an TV-basis mea- 
surement. The X measurement can be effected by first 
applying a Hadamard H operation which rotates from 
the X basis to the Z basis, followed by a computational 
basis measurement. In a ID cluster, measuring the first 
qubit initially in the state teleports the modified 
state X m H ~R z (£)\ip) to its nearest neighbor, where we 

define R a (0 = e~ ?; ^ E with a = x,y,z and £ = TV, Y, Z. 
The measurement outcomes m = and 1 correspond 
to projection into the computational states |0) and |1), 
respectively. Universal single-qubit operations can be de- 
composed into three successive rotations around orthogo- 
nal axes (the Euler angles) , accomplished in the one-way 
model by three successive measurements: 

|V/> = X m 'HR g (0 3 )X m 'HR M {e a )X mi HR t (e 1 )\il>) 
= X^Z^X^HR^RMRz^M, (1) 

where 9' 2 = (-l) mi 2 and 6' 3 = (-I) r ™ 2 <9 3 . Note that 
R x {9) = cxp(— i6X). The Hadamard and the byproduct 
unitaries (the X and Z in Eq. |T]) above) arising from 
measurement outcomes m, = 1 can wait to be applied 
only at the end of the gate teleportation. Crucially the 
angles 9[ for subsequent measurements must be adapted 
based on the result of previous measurements for the 
computation to be deterministic, a property known as 
'feed-forward.' Together with the vertical links in the 
computation-specific cluster state, representing CZ gates 
between computational qubits, measurements in 1WCQ 
generate a set of gates that are universal for quantum 
computation. 



B. Single Atom Rotations 

Rotations on atoms are effected either by a radio- 
frequency pulse resonant on the |0) <-> |1) transition or by 
a Raman transition between these levels using two lasers. 
In either case, the electromagnetic field associated with 
the radiation gives rise to the Hamiltonian in the dipole 
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and rotating wave approximations [2c 



Hi = 



(|l)(0|e^ + |0)(l|e-^) 



(2) 



where f2 and tp are the Rabi frequency and laser phase, 
respectively. The time evolution operator U = er %Ht l h 
is then 



u(nt, ip) = 
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where U(Qt,0) = R x (Qt) and U(Q,t,n/2) = R y (Qt). 

Using the identity H = —iR x (ir)R y (Tr/2) we can write 
HR z (0) = R x (0)H = R x (6 + Tr)Ry(n/2) neglecting the 
overall phase. Because an R y (ir/2) is always present in- 
dependent of the angle 0, one can simultaneously apply 
R v (tt/2) with a wide laser beam to all the atoms in the 
computational cluster state before making any measure- 
ments. In this way only R x rotations need to be applied 
to the qubits before measuring in the computational ba- 
sis: the result of the teleportation is the same as if the 
rotation had been applied after measurement of the pre- 
vious qubit. 

The focused rotation laser is unfortunately not as nar- 
row as the lattice spacing, and will therefore irradiate 
atoms near the qubit of interest. The beam intensity has 
a Gaussian distribution / = Iq exp(— 2x 2 /r 2 ) where x is 
the distance from the center and r is the beam radius at 
which I = l^er 2 . The Rabi frequency fl is proportional 
to the electric field and therefore the square root of the 
intensity tt(x) = Qq exp(—x 2 /r 2 ). All neighboring atoms 
are influenced by the field for the same pulse length t, so 
atoms further from the beam center will undergo a rota- 
tion by an exponentially decreasing angle. 



C. Single Atom Measurements 

Measuring the internal state of an atom is done by 
fluorescence, a technique used for quantum jump detec- 
tion 0, HH, . Atomic measurements of this type have 
been experimentally demonstrated with ions (where the 
technique is also known as electron shelving) and work 
with high probability p7T | . In a fluorescence measurement 
a pulse is applied between the computational state |0) 
and an unstable auxiliary state \2) which rapidly decays 
spontaneously back to |0) (the decay pathway 1 2) — > |1) is 
forbidden) . The atom will repeatedly transition between 
|0) and 1 2) producing many fluorescent photons as long 
as the measurement pulse is active, while an atom in |1) 
will not be affected by the measurement beam. A strong 
measurement corresponds to the projection of the atom's 
internal state wavefunction into one of the computational 
basis states with high certainty. 

While the measurement pulse is active, the Hamilto- 
nian for the system in the dipole and rotating wave ap- 
proximations is given by: 



# =^(|0)(2| + |2)(0|). 



where is the |0) <-> |2) Rabi frequency. The time evo- 
lution of the atomic density matrix is described by the 
master equation [23[ 

p = -^[lf J p] + |(2|0)(2|p|2)(0|-|2)(2|p-p|2)(2|) (5) 

where 7 ^> ft is the rapid decay rate from excited level 
1 2). For times long compared to I/7 one can effectively 
eliminate the population in |2). Neglecting terms pro- 
portional to (O/7) 2 one obtains 
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These equations describe a two-level system with a decay 
of the poi coherence a rate ^-t. A successful measure- 
ment requires the vanishing of the off-diagonal elements 
(coherences) of the density matrix, ^-t > 1, which to- 
gether with O/7 -C 1 requires t ^> fi . The probability 
of a measurement occurring is therefore given by 



P 



1 



(9) 



(4) 



which approaches unity exponentially for high intensi- 
ties, strong coupling, rapid decay times, and long mea- 
surements. 

As was the case for rotations, the spatial profile of 
the measuring laser ensures that atoms near the qubit 
of interest will have a reasonably high (though exponen- 
tially decreasing) probability of being measured. The 
consequences would appear to be disastrous for 1WQC, 
because unwanted projections of qubits in the compu- 
tational cluster will completely destroy the state of the 
logical qubit. Furthermore, an inadvertent measurement 
will be nearly impossible to detect. The fluorescence is 
proportional to the local beam intensity, so that atoms 
projected by the low-intensity wings of the Gaussian will 
emit comparatively few photons. The solution to this 
apparent impasse will be discussed in Section [TV] 



III. CLUSTER STATE CREATION 

The first hurdle facing practical 1WQC in an optical 
lattice system is the formation of the computational clus- 
ter state. The usual approach taken is to first create a 
standard cluster state and then to selectively remove un- 
wanted physical qubits through Z measurements. This 
method is not feasible for wide measuring beams, how- 
ever. There will always be some probability of inadver- 
tently (and unknowingly) measuring the states of nearby 
qubits, thereby removing atoms that should form part of 
the cluster. A more suitable construction is to first rotate 
all computational cluster qubits to |+) and all unwanted 
qubits to |0), and only then perform the controlled phase 
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gates between all neighboring qubits. Qubits in |0) will 
be left untouched by the CZ operations and therefore 
will not be connected to neighboring qubits. 

It is possible to specify the states of individual qubits in 
an optical lattice even for wide rotation lasers that irradi- 
ate multiple atoms. One requires only that the center of 
the beam can be positioned with high accuracy. Consider 
an array of N qubits on which one performs N rotation 
pulses, each centered on a different qubit, with angles n 
about the same axis, where n = 1,2,. ..,N. The beams 
have a Gaussian intensity profile and therefore will ro- 

2 / 2 

tate each qubit m by an angle <j) m = 6 n e x ™- n i r , where 
Xmn = \x n — x m \ is the distance from the central qubit 
n. This will lead to a linear system of equations with N 
unknowns (the n ) that yield the desired phases (the 4>n) 
on each qubit. In general, the linear system takes the 
form 



(a) 



AS 



(10) 



where A is of the form 
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In order for a solution to exist one requires 

det(A) =^6 fc e- fe/r V0, (14) 
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FIG. 1: Building a two-qubit cluster state with three qubits. 
The initial state (a) has all three qubits in state |+); after 
rotations qubit 3 is in state |0), so that entanglement leaves 
the desired final state (b). 



shown in Fig. [1] and one wishes to generate entanglement 
between the qubits 1 and 2 but leave qubit 3 disentan- 
gled, FigQJb). First all three qubits are initialized to |+) 
with a very wide beam (which assumes negligible spatial 
variation in the intensity). Three R v {9i) pulses are then 
applied, each one centered on a different qubit, such that 
qubits 1 and 2 remain in |+) while qubit 3 returns to |0). 
The linear system is: 



l:6 l 

2 : e- 1 / r ' 2 9 l 

3 : e~ 1 / r2 9 l 

which has the solution 
0i = 



+ e 



-2/r 2 



e- 2 ^02 



7T e 



-1/r 2 



2 1 - e- 2 A 2 ' 
0; 

7T 1 

~2 1 -e- 2 /- 2 



0; (15) 
0; (16) 
-tt/2, (17) 



(18) 
(19) 
(20) 



In this particular case only two rotations are actually 
needed. Note that the solution depends on the beam ra- 
dius r only in terms of an (experimentally) adjustable 
parameter, as long as r remains finite. For a larger num- 
ber of qubits and a wide beam, each rotation will affect 
many qubits across the lattice, but as long as there are as 
many distinct operations as qubits, each applied rotation 
may be chosen such that the overall phase on each qubit 
is exactly as desired. 



where both the k (sums of distances squared in the x 
and y directions of the cluster) and the bk (combina- 
toric coefficients) arc integers. Suppose one assumes that 

det(A) = 0. This is equivalent to stating that e _1 ' r is 
a root of a polynomial with integer coefficients. For ra- 
tional — 1/r 2 , however, e -1 / r is transcendental which by 
definition cannot be the root of a polynomial with integer 
coefficients 28]. It follows that det(A) is non-zero and a 
solution exists for the linear equations (|10[) . Therefore, 
with one beam centered on each of the ./V atoms it is 
possible to precisely control the rotation applied to each 
individual atom. 

Consider a simple example in order to demonstrate the 
procedure. Suppose there are three qubits oriented as 



IV. ID CLUSTER COMPUTATION 
A. Offset Measurements 

The simplest way to reduce the probability of inadver- 
tently measuring qubits in an entangled chain is to offset 
the center of the measuring beam. Consider the very 
first qubit in a one-dimensional cluster, Fig. [5] To the 
left of this qubit there is empty space and to the right 
are connected qubits. We wish to measure the first qubit 
with high probability while preserving the states of its 
neighbors. This can be accomplished by taking advan- 
tage of the Gaussian shape of the laser pulse: shifting the 
beam to the left while increasing its maximum intensity 



5 




FIG. 2: The advantage of offsetting Gaussian pulses. The 
centered beam M c has the same intensity on the target qubit 
as a higher-intensity beam M whose center is displaced from 
the target qubit by a single lattice spacing to the left, but M 
has an exponentially reduced intensity on the first neighbor 
to the right. 



reduces its amplitude on the neighboring qubits without 
altering its effect on the target atom. 

Suppose that one requires intensity I = It on the 
first (target) qubit. Shifting the beam center n lat- 
tice spacings a to the left yields intensity on the tar- 
get I = I n e" 2 '™"^ / r , so that the intensity maximum 
of the beam must increase exponentially in the offset: 
I n = I t e 2( < na } / r . The intensity on the first neighbor 
to the right of the target If nevertheless decreases expo- 
nentially in the offset: if = I^i™) 2 /r 2 e -2(n+l) 2 a 2 /r 2 = 
j te -2(2n+i)a /r _ jf ^he desired probabilities of measur- 
ing the target and nearest-neighbor qubits are p t and pf , 
respectively, then from Eq. ((HJ) one obtains the required 
offset 

(21) 



n = — - In 

4a 2 



In(l-Pt) 



ln(l - p { ) 



For example, if p t = 0.99 andpf = 0.01 were desired, then 
with r = 4a one would require n ps 24. With r = 10a 
the offset jumps to approximately 153a. The probability 
of simultaneously measuring the second nearest neighbor 
p s is truly negligible: 



p B = 1 - exp 



2 /r 2ln 2 (l — pt) 



-4a 2 /r 



In(l-pt) 



(22) 



For the same choices of pt and pf above, one obtains 
p s a 10" 5 . 

The discussion above has assumed that the target 
qubit is located at the left-most edge of the chain, with 
its neighbors to the right. In the most general case of 
an offset measurement, however, a significant number of 
qubits M to the left of the target will also be measured 
at each stage of the computation (M ps 2n where na is 
the offset distance discussed above) . Because all of these 
were target qubits during previous steps in the one-way 
computation, their subsequent projections into computa- 
tional basis states during the measurement of the target 
qubit will not affect the logical states of any of the cluster 
qubits. That said, rotation pulses will have been applied 
to these qubits by the wide lasers during the computa- 
tion, so their measurement outcomes are not predeter- 
mined. 



One must be able to clearly distinguish the fluores- 
cence of the target atom from that of all the qubits to 
its left when performing a measurement with an offset 
beam. The simplest solution is to re-set the states of 
all M non-cluster qubits that will be re-measured to the 
non- fluorescing state |1). This can be accomplished with 
the same method as described in Sec. Mil to initialize all 
of the qubits' states. The phases ipi accumulated on the 
M relevant qubits by previous rotations are known, so 
that one can choose <f>i = — tpi + tt as input to the linear 
system of equations to find the 0j. Note that the appli- 
cation of the M rotations will also build up an undesired 
phase on the target qubit, so this qubit needs to be in- 
cluded in the linear system, with ^target = 0. Thus, M+l 
operations are required for each target measurement. In 
addition, the phases of 2n qubits need to be calculated 
classically at each stage; together with the time for the 
classical solution to the 2n + 1-dimensional linear system, 
each measurement has a classical overhead that scales like 
4n 2 . 



B. Measurement Protocol 

As discussed in detail above, with wide beams one is li- 
able to make inadvertent measurements of atoms close to 
the target qubit, though the probability of this occurring 
can be made arbitrarily small in principle by offsetting 
the center of the measuring beam from the target. If only 
one of the measured qubits to the right of the target is 
in a state outside the XY plane prior to its measure- 
ment, the logical qubit will undergo a non-unitary (and 
therefore uncorrectable) transformation. The most ex- 
treme case is a Z measurement, which unentangles the 
measured qubit from the rest of the cluster and destroys 
the entanglement resource needed for quantum telepor- 
tation. Even if all measured qubits' states were initially 
in the XY plane, the absence of feed-forward makes it 
unlikely that the desired unitary transformation will be 
implemented. 

There are two cases when feed-forward is not required, 
however. The first is if all of the measurement out- 
comes are zero, in which case there are no byproduct 
operators to commute through the applied rotations, and 
therefore no compensation of the choice of measurement 
angle. Unfortunately, the likelihood of this occurrence 
decreases exponentially with the number of measured 
qubits, and in any case the results of inadvertent mea- 
surements arc by definition unknown. The second case 
corresponds to teleportation of single qubit Clifford oper- 
ations, which byproduct operators can commute through 
without changing the operation [29l |. It is impossible, 
however, to effect universal quantum computation using 
Clifford operations alone. 

Below we describe in detail a protocol that elimi- 
nates errors arising from inadvertent measurements of 
near neighbors, with linear overhead in terms of ancillary 
physical qubits and operations. Feed-forward is normally 
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accounted for by applying a corrective rotation depend- 
ing on the result of the previous measurement. As dis- 
cussed in Sec. Ill A| this corrective rotation need only be 
applied when that result is one. When mi = the gate 
teleportation is successful even if an inadvertent measure- 
ment of a neighbor has occurred. When mi = 1 there will 
be some probability p that an error has occurred, but the 
negative consequences of this can be avoided by imple- 
menting feed-forward in an alternative manner. Rather 
than applying a corrective rotation on the next measure- 
ment, we resign ourselves to rotating in the wrong direc- 
tion. Instead, the sequence of measurements is extended 
by a number of qubits on which the phase error can hope- 
fully be corrected. We first describe in detail the protocol 
to compensate for inadvertent measurements of the qubit 
that is the nearest-neighbor of the target; the extension 
to the case where an arbitrary number of qubits is inad- 
vertently measured is discussed at the end. 

The goal is to implement three successive Euler rota- 
tions HR z (ai), with the ai unequal angles. As discussed 
in Sec. Ill Al these constitute a universal singlc-qubit uni- 
tary on the computational qubits. After entanglement, 
R v {-k/2) is applied to all of the qubits so that future rota- 
tions can be performed only about the X axis. To effect 
the first two Euler rotations, one applies R x rotations 
such that the first and second qubits are chosen to have 
a total angle of ai + tt and a 2 + n, respectively. This can 
be accomplished by applying two rotation pulses in anal- 
ogy to the cluster-carving protocol described in Sec. IIIII 
except now only with R x rather than R y rotations. The 
first pulse rotates the first qubit by an angle 9\ and the 
second pulse rotates the second qubit by 6* 2 , where the 6i 
can be chosen by solving the linear system 



a 2 +TT = e~ 1/r2 9! 



(23) 



(24) 



with r the width of the rotation beam. Solving for 9\ 
and 02 one obtains 



*i + 7r - e 1 / 1 ' 2 [a.2 + 7r) 
1 - e" 2 /- 2 ' 



6i 



a 2 



+ TT-e- 1 / r2 (a x +tt) 



1 - e- 2 / r ' 2 



(25) 



(26) 



These pulses will apply undesired unitarics to the other 
qubits in the chain, but these can be compensated for 
later. The total applied rotations are then 



R x {a x +ir)Ry{-K/2) 
R x (a 2 +Tr)R v (Tr/2) 

R x (e- i ' r2 e l + e- 1 / r2 e 2 )Ry(TT/2) 



i.e. total unitaries of HR z {a\) and HR z (a2) have now 
been applied to qubits 1 and 2, respectively. 

Measurement of the qubit 1 teleports the logical state 
X' mi HR z (ai)\ip) to the second qubit. Taking into ac- 
count the finite probability of an inadvertent measure- 
ment of the second qubit, the logical state (defined as the 
state of the qubit immediately to the right of the atom 
just measured, not including as-yet unmeasured qubits 
further to the right or single-qubit unitaries that have 
resulted from previous rotations) is the density matrix: 

Pl = (\-p)\\X m -HR z { ai )m\ 

+ p\ \X m * HR Z (a 2 ) A™ 1 HR z ( ai )\m, 

= (l-p)\\X^HR z ( ai )\m 

+ p\\X m -Z m 'R x [(-l) m 'a 2 }R z ( ai M\\, (27) 

where HIV') 1 1 = IV'XV'I; P = Pu an d the identities 
R z (a)X = XR z (—a) and R x (a) = HR z (a)H have been 
used to derive the final expression. Note that the fidelity 
can never be exactly unity during the computation be- 
cause there is always the probability p of inadvertently 
measuring the adjacent qubit. 

Now the second qubit must be measured. If the second 
qubit had already been inadvertently projected, then the 
outcome of the second measurement will definitely be m 2 ; 
otherwise the value of 11J2 will be determined presently. 
If toi = 0, one could carry out this measurement by leav- 
ing the first qubit in its |0) state and distinguishing the 
second qubit 's measurement outcome by an increased or 
unchanged level of fluorescence compared to measuring 
the first qubit alone. In keeping with the approach out- 
lined in the above section, however, prior to measurement 
one should rather apply a total rotation of 7r to the first 
qubit to rotate its state into |1) (which produces no fluo- 
rescence signal), and a rotation of to the second qubit 
to preserve its state. 

Neglecting for the moment the possibility of inadver- 
tently measuring the third qubit upon measurement of 
the second qubit, the result would be the pure logical 
state X m2 Z mi R x [(-l) mi }R z ( ai )\ij). If mi =0, then the 
second Euler rotation will have been correctly effected, 
and one could simply apply an angle a 3 + n to qubit 3 
to effect the last Euler rotation, completing the general 
singlc-qubit unitary. (Recall that with the wide beam 
one also needs to ensure that all previously intention- 
ally measured qubits are rotated into |1), and that the 
accumulated phases on atoms beyond qubit 4 are noted 
for future reference). Now explicitly including the pos- 
sibility of inadvertently measuring the third qubit, after 
measurement of qubit 2 the logical state would be 

P 2 = (l-jOII* ma ik(«a)fl*(«iMII 

+ pWX^HR^X^R^R^mW 

= (l-rfllA^M^Ml^H 

+ P \\X m °Z m *HR z {(-l) m *a 3 ]R x (a 2 ) 

xR z ( ai )\iP}\\. (28) 

Again, if m 2 = 0, then the third Euler angle will have 
been applied correctly, and measurement of qubit 3 will 
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complete the unitary. The rotation on the fourth qubit 
is arbitrary, so let's set it to tt. The logical state is finally 

p 3 = (l-p)\\X m3 HR z (a 3 )R x (a 2 )R z ( ai M\\ 

+ p\\X mi HX ms HR z (a 3 )R x (a 2 )R z (a 1 )\tP)\\ 

= (l-p)\\X^HR z (a 3 )R x (a 2 )R z ( ai m\\ 

+ p\\X m *Z m 3R z (a 3 )R x (a 2 )R z ( ai )\i;)\l (29) 

It is important to note that the overall fidelity has not 
decreased from the state (|27|) after the first measure- 
ment. Each successive measurement effectively purifies 
the logical state into the same sum of two density ma- 
trices. The right-most cluster qubit is well-defined and 
so the very last measurement can be chosen to be un- 
ambiguous, yielding unit fidelity at the end of the full 
1WQC. 

The scenario above assumes that m, = for all i, but 
the probability of this occurring decreases exponentially 
in the number of measured qubits. Thankfully, post- 
selecting on these outcomes is not necessary with our 
protocol. Consider the case mi = 1. Measuring qubit 
2 will result in the incorrect Euler angle — a 2 . This sce- 
nario is quite different from standard 1WQC, whereupon 
measuring mi = 1 the opposite angle would be chosen for 
the measurement of qubit 2 to immediately correct the 
feed- forward error. In the present case qubit 2 may have 
already been inadvertently measured, and it is too late to 
correct the error in this manner. Rather, measurements 
are performed on the next pair of qubits in an attempt 
to rotate the Euler angle back to the desired value. Ro- 
tations of tt and —2a 2 + tt are applied to qubits 3 and 
4, respectively. Measuring qubits 2 and 3 then yields the 
logical state 

p 3 = (l-p)\\X m3 HX m2 ZR x (~a 2 )R z ( ai M\\ 
+ p\\X" u HR z (-2a 2 )X ms HX m2 ZR x {-a 2 ) 

xiZ,(ai)fo&>|| 
= (1 -pjWX^Z^-XHRzi-a^R^a^W 
+ p\\X mi Z m3 X m2 ZR x [{-l) m3 2a 2 ]R x {-a 2 ) 

x^(ai)|V>||. (30) 

If m 3 = 0, then the logical state becomes 

P3 = (l-p^Z^-XHR^-a^R^a^M 

+ p\\X m ^-ZR x (a 2 )R z ( ai m\l (31) 

so that the correct second Euler would be applied upon 
measurement of qubit 4. It suffices to rotate qubit 5 
into a 3 + tt in advance, and continue as if no error had 
occurred. If rather m 3 = 1, then the applied Euler an- 
gle would be —3a 2 , and another qubit pair need to be 
inserted (the second with angle 4a 2 + tt), etc. The prob- 
ability of not obtaining a outcome on the relevant odd- 
numbered qubit (and therefore not applying the correct 
Euler rotation) decreases exponentially in the number of 
attempts. 

Following this example above, if now m 3 = one can 
choose qubit 5 to have the phase (— l) Tn2 a3 + tt, so that 



measurement of qubit 4 yields the logical state 

p 4 = {l-p)\\X mM ZR x (a 2 )R z ( ai m\\ 
+ p\\X m5 HR z [{-l) m -a 3 ]X mM ZR x (a 2 ) 

xR z ( ai )\M 
= (1 -p)\\X m ^ m -ZR x (a 2 )R z (a x )W)\\ 
+ p\\X m "Z mM XHR z [(-iy^a 3 ]R x (a 2 ) 

xR z ( ai )\M- (32) 

If 777,4 = 0, measurement of qubit 5 yields the desired 
unitary; otherwise, additional pairs of qubits need to be 
appended to qubit 5 to rotate the Euler angle back to 
the desired value, as described above. 

We now summarize the protocol. Each HR z (ai) oper- 
ation in principle requires only one measurement on each 
physical qubit, all prefaced by the appropriate rotations 
on nearby qubits. However, if the sum of certain previous 
measurement outcomes is odd, one needs to correct for 
errors due to possible inadvertent measurements (even if 
they haven't occurred!) by inserting additional pairs of 
measured qubits. Four qubits is the minimum number 
for a universal singlc-qubit unitary, though in practice 
the number could be considerably higher, depending on 
the measurement outcomes. 

While this protocol protects against inadvertent mea- 
surements on the qubit immediately to the right of the 
target, accidental measurement of the next-nearest neigh- 
bor could be catastrophic for the overall 1WQC. With 
offset measurements the probability of this occurring is 
small, but for a long computation the number of such 
events can become appreciable. A protocol protect- 
ing against inadvertent measurements of many neighbors 
avoids the requirement of very large laser offsets with 
wide beams. It could also provide a way to protect quan- 
tum information against re-absorption by distant atoms 
of scattered photons emanating from the fluorescing tar- 
get. 

To protect against more long-range inadvertent mea- 
surements requires additional 'buffer' qubits rotated by tt 
that effect only Clifford gates. Consider the case where 
one intends to compensate for inadvertently measuring 
two neighbors. As above, all qubits initially undergo 
R y (rr /2) rotations. The first two qubits undergo R x rota- 
tions by a±+TT, a 2 + n as before. Now the possibility that 
the third qubit might be inadvertently projected sug- 
gests that it be rotated by tt prior to any measurements. 
Again, the outcome of the first measurement m\ deter- 
mines if the second Euler rotation will be implemented 
correctly. If mi = 0, then one can proceed directly with 
implementing the third Euler rotation. One cannot sim- 
ply rotate qubits 4 and 5 respectively by a 3 + tt and tt, 
however. The combination of tt and a 3 + tt on qubits 3 
and 4 eliminates the Hadamard gate needed for the next 
Euler rotation to be around the z axis. One therefore 
needs qubits 4 through 6 rotated by tt, (— l)™ 2 a3 + tt, 
and tt, respectively. Again, 7714 determines the success of 
the last Euler rotation. If 7714 = then qubit 7 needs to be 
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set to 7r and the fifth qubit measured to effect the single- 
qubit unitary. Thus, a minimum of seven qubits are re- 
quired to protect against two-qubit inadvertent measure- 
ments. If either mi = I or 7714 = 1, additional qubit pairs 
must be added until the Euler angle is rotated back to 
the desired value, just as in the case considered in detail 
above; in the former case, for example, qubits 4 and 5 
would be rotated by — 1012 + t and tt, respectively. 

It should now be apparent that to protect against in- 
advertent measurements on m qubits to the right of the 
target requires on the order of 3m physical qubits in order 
to effect a universal singlc-qubit gate on computational 
qubits. Because feed- forward errors occur with probabil- 
ity 1/2, the number of qubits needed to implement the 
correct Euler angle is of order m. A much larger number 
of measurements may be required in practice to correct 
an Euler angle error, but the probability of continued 
failures becomes exponentially smaller in the number of 
measurements. 

Compensating for inadvertent measurements on mul- 
tiple qubits is not equivalent to making actual measure- 
ments on multiple qubits. The former assumes that flu- 
orescence from inadvertently projected qubits is not ob- 
servable while the latter assumes that it is. To simulta- 
neously measure multiple qubits one could in principle 
make use of the Gaussian profile of the measuring pulse: 
the qubit fluorescence signals are proportional to the lo- 
cal pulse strength. The ability to clearly identify the 
signal strength with an atom's position degrades rapidly 
with the pulse width, however, particularly considering 
the inherent noise associated with photon counting. 



V. 2D CLUSTER 
A. Scheme 

The full two-dimensional protocol works much the 
same way as the ID case discussed in Sec. IIVI In the 
simplest scheme, the beams arc offset horizontally as be- 
fore, but now are centered vertically. The cluster must 
be initially carved so that pulses measuring the physical 
qubits on one chain have a negligible probability of pro- 
jecting the states of physical qubits on an adjacent chain. 
The horizontal chains representing computational qubits 
must be sufficiently well-separated vertically. Naively one 
might assume that the chains would need to be separated 
by a distance at least greater than na, with n the offset 
defined in Eq. (f2"Tj) . 

In practice the chains can be positioned much closer 
together, because of the circular spatial 'footprint' of the 
measuring pulse. Consider a beam that is centered ver- 
tically on a given chain but is offset horizontally by na. 
If the intensity on the first qubit is Iq, then the first 
qubits on chains above and below will experience inten- 
sity 7 e- 2 (" la ) 2 / r2 , where ma is the separation between 
horizontal chains. If the desired probability of projecting 
the first qubit on an adjacent chain is p m , then the chain 



separation is 




In(l-jJt) 



In (I - Pm) 



(33) 



Thus the inter-chain separation grows more slowly (by 
a factor of r / a) than the offset required to minimize in- 
advertent measurements on adjacent qubits on the same 
chain, Eq. (f2"Tj) . If for example one chooses p t = 0.99 and 
a tiny value p m = 1 CP 10 , then m fw 3.5(r/a); with r = 4a 
and r = 10a one obtains m ~ 14 and m ~ 35, respec- 
tively. This separation is adequate for all measurements 
on all chains (momentarily postponing issues related to 
links between horizontal chains, discussed in detail be- 
low). One simply requires that all the measured qubits 
in the various chains have the same horizontal coordinate 
(column index). 

Recall that prior to each measurement in the ID pro- 
tocol, one needs to perform of order 2n rotations. These 
prepare the target qubit and those to its right that might 
be inadvertently measured, and restore qubits previously 
measured and soon-to-be- remeasured to the state |1). In 
2D, the number of measured qubits, and therefore the 
number of preparatory rotations, at each stage is of or- 
der irn 2 . This requires a solution to a n 2 x n 2 linear 
system of equations whose cost scales as n 4 , a signifi- 
cant classical overhead. Note that the beam irradiating 
a given chain impinges on chains above and below. So 
the necessity of performing preparatory rotations implies 
that measurements of qubits on chains separated by dis- 
tances smaller than na cannot be done in parallel, unless 
the solution of the relevant linear equations involves all 
rows of physical qubits. In practice, this apparent lack of 
operational parallelism is no serious impediment: in the 
worst-case scenario considered above with r = 10a, there 
are only 4.4 chains per full offset. 

The main extension of the ID protocol to the 2D clus- 
ter is how to implement the teleportation primitive on 
the links between horizontal chains, in order to simulate 
an entangling gate between computational qubits. As 
was discussed in Sec. Ill Al no feed-forward is required 
when single-qubit measurements teleport Clifford gates. 
Because the CZ gate belongs to the Clifford group, with 
suitable preparation of the relevant qubits (i.e. applica- 
tion of total rotation of tt) the inter-chain links (hith- 
erto referred to as simply links) should be unaffected by 
inadvertent measurements. That said, the outcomes of 
measurements on links need to be unequivocally known, 
since they will affect the feed-forward criteria for future 
measurements on chains. 

Making well-defined measurements on qubits in verti- 
cally oriented links is extremely difficult, however. Con- 
sider the situation where the target qubit on a given chain 
is the first qubit of a vertical link. Using Eq. ([5]), the 
probability of measuring the first link qubit is close to 
p t . Indeed, Eq. (|3"3"|) states that m qubits along the link 
from the target will be measured with probability p m 
or greater. If we wish that this probability be at least 
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smaller than the probability of inadvertently measuring 
a qubit to the right of the target along the chain p m > < pf , 
then the number of link qubits that will be measured is 
m' = y/2n + 1 « m/2 using Eq. (|2"Tj) . Thus, the states 
of individual link qubits will be impossible to uniquely 
determine. 

The solution to this apparent conundrum is to form 
inter-chain links that are essentially diagonal, with the 
'zig-zag' pattern shown in Fig. [3J The first two qubits in 
the link must be vertically oriented. Recall that the clus- 
ter is formed by first initializing all cluster qubits to |+) 
followed by the CZ gate, and one needs to avoid forming 
a box graph. With this geometry, the maximum number 
of inadvertent measurements will be the number of qubits 
in the column to the right of the target. It is clear from 
the figure that this number will often be three: one on 
the chain and two nearby in the link. More importantly, 
the maximum number of qubits intentionally measured 
simultaneously will also be three. Even without the abil- 
ity to spatially resolve the signal emanating from three 
atoms, however, the individual states of the three qubits 
can still be obtained with minimal operational overhead 
as shown explicitly below. 



B. Protocol 

The 2D protocol differs primarily in that simultane- 
ous measurements of chain and link qubits must be per- 
formed. Suppose the first target is in column 1, on the 
upper chain just left of the first link qubit, as illustrated 
inFig.[3l It is assumed to be in the state = a\0)+[3\l) 
with a, (3 complex coefficients satisfying \a\ 2 + \f3\ 2 = 1. 
If all entangled qubits are first rotated by n, after mea- 
surement of qubit 1 the resulting logical state is 



Pi 



« {l-pt-f4-pt)\\X m *H\ 

+ p 2 t \\X^H 2 X^H x H 2 X^H r \ 

+ pfWX^HiX^HsX^HiHzX^H^W 
= (l-p f -p?-p f 3 )||A mi F|V>>|| 

+ Pf \\x^z^HS)\\ + P i\\xpx™°zrm\ 

+ p!\\XpZ^X^Z^H 2 \i>)\\- (34) 



This expression assumes that inadvertent measurements 
are possible only on qubits in column 2 (labeled qubits 
2, 3, and 4 in Fig. [3]), in the vicinity of the upper chain. 
The small differences in the probabilities of inadvertently 
measuring qubits down the link are neglected for clarity, 
and are all set to pf. Note that to define the current 
logical state one can ignore measurements beyond the 
qubit that connects to the original logical state. Also, 
in practice one can ignore terms with prefactors higher 
than p\ , since these are much smaller than the probabili- 
ties of inadvertently measuring second-nearest neighbors 
which arc ignored in the current treatment. In any case, 
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FIG. 3: The full 2D protocol. Initially, qubit 1 in the first 
column of the upper chain is measured by a beam that is offset 
horizontally but centered vertically (the right-most edge of the 
beam is labeled mi). Immediately thereafter, the first qubit 
on the second chain is measured with the same horizontal 
offset, as is the first in the third chain, etc. Next, qubits 2, 
3, and 4 in the second column are simultaneously measured 
(7T12), as is the second qubit in the lower chain, and so on. 
Measurements continue in this fashion until the logical state 
is entangled between the two chains. 



the outcomes of these inadvertent measurements will be 
determined shortly by intentional measurements. 

The most important new feature in the density matrix 
(|3~4")l is the appearance of the state \ip) = a\00) + /3|11). 
When a measurement is performed on the chain qubit 
that is also the first link qubit, the chain and link 
qubits become inherently entangled. The logical state 
must therefore be written in a two-particle basis with 
|0) i— ► |00) and |1) 1— > |11), with the first and second 
elements representing the chain and link qubits, respec- 
tively. Of course, this is how the link is able to simulate a 
CZ gate between two chains: upon measurement of the 
last link qubit, which is located on the second chain, the 
two qubit basis corresponds to chain 1 and 2. The sub- 
scripts on the single-particle operators in Eq. (j34|) there- 
fore correspond to the qubit index in the two-particle 
basis. 

After rotating all nearby unmeasured cluster qubits 
by 7r, the measurement of atoms in the second column 
(qubits 2-4 in Fig. [3]) can be performed. The outcome is 



p 2 a (l-2p { -2p 2 )\\X™*Z™ 3 X™Z^H 2 



+ ptWX^X^Z^Z™ 2 X™ 1 H 2 Hi 
+ p f ||A 2 m6 Z 2 m4 A™ 3 A™ 2 Z;' 



p\ 1 | A™ 7 Z™ 6 A 2 " 14 Z™ 3 X™ 2 Z\ ni H 2 \i>}\\ 
p 2 1 1 A™ 6 A™ 5 Z™ 4 A™ 3 Z™ 2 A™ 1 Hx\i>}\\, 



(35) 



10 



ignoring terms of order pj? and higher. How is it possible 
to uniquely determine the outcome \m2m3m4)? If no sig- 
nal is obtained, the three-qubit state is definitely |111). 
If a weak signal indicates only one atom is in |0), then 
the possible outcomes are |011), 1 1 1 ) , and 1 1 10) . Subse- 
quent rotation of the qubit 2 by tt will either yield 1 1 1 1} 
or two qubits in |0), which can be easily verified by re- 
measurement of the three qubits. If the latter outcome is 
obtained, then re-rotating qubits 2 and 3 both by tt will 
unequivocally determine the correct initial output. The 
same approach works for an original measurement indi- 
cating two |0) outcomes, and trivially for |000) or |111). 
Even if the relative signal strengths corresponding to dif- 
ferent numbers of |0) outcomes cannot be distinguished, 
there are only seven different combinations that need to 
be attempted (by repeated rotations and measurements) 
in order to convert the output to |111) and thereby elim- 
inate the fluorescence. 

Measurements continue in this columnwise fashion un- 
til the link connects with the lower chain. The signal 
strength from link qubits on measurements of chain 1 
qubits will gradually decrease, becoming undetectable 
soon after the midpoint between chains. Before this col- 
umn is reached, measurement of chain 2 qubits will pick 
up signal from the link qubits, maintaining the flow of 
information down the link. If the fluorescence noise is 
too great, then one would simply move the center of the 
measuring beam vertically between the two chains to im- 
prove the signal. 



VI. SUMMARY AND DISCUSSION 

We have described in detail a protocol for deterministic 
quantum computation in the one-way model, where mul- 
tiple physical qubits are simultaneously measured. This 
approach is well-suited to implementation with ultracold 
atoms confined in optical lattices, where the measuring 
beam is much wider than the separation between qubits. 
A central assumption of the approach is that the cen- 
ter of the lasers can nevertheless be positioned with high 
accuracy. Our strategy makes use of the Gaussian pro- 
file of lasers used in these systems for performing both 
rotations and measurements. For rotations, the Gaus- 
sian profile allows for a unique rotation on each clus- 
ter qubit to be performed by superimposing rotations 
on all physical qubits. This enables the direct genera- 
tion of the desired computational cluster state, and to 
telcport the desired unitaries during the one-way com- 
putation. For measurements, the laser center is chosen 
to be far removed from the qubit of interest in order to 
minimize the possibility of making inadvertent measure- 
ments on nearby qubits. With additional overhead in 
terms of physical qubits, however, the protocol can ac- 
commodate an arbitrary number of unknown outcomes 
of measurements accidentally performed on nearby con- 
tiguous qubits. 

Each measurement, which effects a gate on the com- 



putational qubits, is performed in four steps. First, all 
qubits that will be irradiated by the measuring beam, as 
well as several qubits beyond, must first have their phases 
suitably prepared. This is accomplished by applying a ro- 
tation pulse centered on each physical qubit. Second, the 
relevant cluster qubits are measured. Third, if more than 
one qubit is intentionally measured, then these need to be 
further rotated and subsequently re-measured. Fourth, if 
previous measurement outcomes will incorrectly apply a 
future rotation due to an inadvertent measurement, then 
this error needs to be fixed by a finite set of future mea- 
surements. 

Because telcportation errors caused by inadvertent 
measurements need to be fixed by extending the num- 
ber of measurements of qubits in the horizontal chains, 
it is not known a priori how many chain qubits will be 
needed to effect the desired unitary. This is somewhat 
problematic because the computational cluster, i.e. the 
entire structure of horizontal chains and vertical links, 
needs to be formed in advance of any measurements. The 
simplest solution is to make all the inter-link distances 
sufficiently long that the total probability of performing 
the correct unitary will be high. The probability of con- 
tinued failure after £ attempted corrections is l/2 e+1 . If 
the desired probability of success is 1 — e then the num- 
ber of correction attempts must be I > log 2 (l/e) — 1. 
In order to protect against inadvertent measurements of 
m neighbors, an arbitrary unitary will require a chain of 
3(lm + 1) qubits for even m or 3[l(m + 1) + 1] for odd m. 
Once the correct unitary is accomplished, measuring the 
remaining qubits will teleport Clifford gates with trivial 
consequences. 

An attractive feature of this scheme is that in princi- 
ple it could partially mitigate the second most significant 
impediment in the implementation of 1WQC with ultra- 
cold atoms in optical lattices: stray fluorescence from 
measurements accidentally projecting distant qubits. Be- 
cause photons are emitted in all directions, the probabil- 
ity of inadvertently measuring qubits a distance d from 
the target decreases like d~ 2 . Simply buffering the chains 
with 32 extra qubits (each rotated by tt to teleport Clif- 
ford gates) reduces the deleterious effects of stray light 
by a factor over 1000. 

While the 1WQC scheme outlined in this manuscript 
can account for the possibility of having inadvertently 
measured qubits, it is currently not designed to protect 
against the small but finite possibility that the target is 
not actually projected by the measurement. A central 
assumption throughout has been that the laser coupling, 
measurement time, and the decay rate are all sufficiently 
large that an intentional measurement is performed with 
high certainty. The issue might seem irrelevant because 
any given target is consistently re-measured as the pro- 
tocol moves down a chain or inter-chain link, so it will 
quickly become projected if it was not initially, and any 
error in assumptions will be made manifest. The prob- 
lem is that rotating a qubit that is believed to have 
been already projected, but is in fact entangled with 
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yet-unmeasured cluster qubits, will possibly tclcport an 
unknown Pauli byproduct gate on subsequent measure- 
ments. This issue is however beyond the scope of the 
current calculations, and will be the focus of future work. 
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